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ABSTRACT:  The  aim  of  this  paper  is  to  obtain  accurate  reference  data  for  the  relatively  large  and 
realistic  reflector  antenna  systems.  Previously  it  has  been  done  for  a parabolic  reflector  antenna  in  E-case 
and  now  it  is  similarly  performed  for  the  H-polarization  case.  Directive  primary  feed  is  modeled  by  the 
complex  source  point  method  and  the  relative  accuracy  of  the  results  is  verified. 

INTRODUCTION:  Reflector  antennas  play  an  important  role  in  the  modern  communication  systems, 
and  hence  attract  research  efforts  for  efficient  numerical  analysis  and  optimization.  Although  two- 
dimensional  (2D)  modeling  relates  only  to  cylindrical  reflectors,  it  is  still  interesting  in  terms  of  advanced 
solutions  and  some  practical  applications.  The  scattering  from  a 2D  reflector  can  be  very  accurately 
simulated  by  analytical-numerical  techniques.  One  of  them  is  the  Riemann-Hilbert  Problem  (RHP) 
technique,  which  reduces  the  canonical  circularly  curved  screen  problem  to  a Fredholm  2nd  kind  matrix 
equation  after  explicit  inversion  of  the  static  part  [1],  In  [2],  we  had  combined  this  solution  with  the 
complex  source-point  (CSP)  method  that  enabled  us  to  model  realistic  feed  patterns  in  the  both 
polarizations.  However  most  of  the  realistic  reflector  surfaces  are  parabolic  ones.  This  necessitates  a 
modification  of  RHP  technique  to  arbitrarily  curved  strips.  In  [3]  we  performed  the  analysis  of  2D 
reflectors  of  conical-section  profiles  in  the  E-polarization  case.  Here  we  used  a dynamic  singularity 
extraction  and  inversion,  as  proposed  in  [4],  to  convert  a logarithmic-singularity  integral  equation  to  the 
Fredholm  2-nd  kind  matrix  equation.  Furthermore  the  application  of  the  FFT  algorithms  enabled  us  to 
treat  electrically  large  geometries.  In  the  H-polarization  case,  the  basic  integral  equation  has  a strong 
singularity.  It  is  worth  to  mention  that  similar  problem  was  considered  in  [5].  Here,  the  integral  equation 
was  reformulated  into  a dual  series  equation  (DSE)  with  trigonometric  kernel.  Then  it  was  regularized  by 
using  a specialized  method  developed  earlier  for  the  equations  with  Jacobi  polynomials.  Later  on,  in  [6]  it 
was  proposed  to  regularize  the  same  DSE  by  the  RHP  technique.  In  the  both  studies,  although  a method 
was  outlined,  no  numerical  data  were  given  for  arbitrary  profiles. 

In  this  study,  we  regularize  Electric  Field  Integral  Equation  (EFIE)  obtained  for  the  parabolic  reflector 
antenna  in  the  H-polarization  case.  EFIE  is  discretized  by  transforming  it  into  the  spectral  domain  and 
also  written  as  DSE.  To  make  regularization  the  static  part  is  inverted  analytically  by  the  RHP  technique 
since  the  static  part  of  DSE  constitutes  a proper  canonical  form.  Furthermore  we  can  say  that  this  is 
equivalent  to  the  semi-inversion  of  the  original  IE  containing  hyper-type  singularity  in  the  kernel  static 
part.  In  the  details  we  used  the  FFT  algorithm  to  perform  computations  of  the  double  Fourier  Series  (FS) 
coefficients  more  rapidly.  All  these  provide  us  to  solve  relatively  large  reflector  antennas  compared  to  the 
previous  results  with  the  semi-inversion  techniques.  This  is  the  continuation  of  the  work  in  [2,3]  and  the 
results  can  be  used  to  check  the  accuracy  of  the  asymptotic  and  numerical  techniques. 

FORMULATION:  The  geometry  of  the  problem  can  be  defined  so  that  infinitely  thin  PEC  curved 
reflector  surface  having  parabolic  profile  is  illuminated  by  a directive  feed  (Figure  1).  Parabola  is  a 
special  case  of  an  arbitrary  profile.  In  the  real  space,  the  feed  is  located  at  the  geometrical  focus  of  the 
parabola.  To  perform  regularization,  parabolic  reflector  surface  is  completed  to  the  closed  contour  C by  a 
circle  having  its  origin  on  the  x-axis.  Its  radius,  a,  is  chosen  so  that,  at  the  connection  points,  the  slopes  of 
the  parabola  and  circle  are  matched.  So  the  contour  first  derivatives  are  continuous  and  the  discontinuities 
in  the  second  derivatives  are  finite.  The  boundary  value  problem  can  be  stated  as  the  satisfaction  of  the 
Helmholtz  equation,  Sommerfeld  radiation  condition,  and  PEC  boundary  condition  on  reflector  surface 
M.  By  usung  the  ftee-space  Green’s  functions  in  2D,  i.e.  GA 0 =G'!>„=i/4  Hj!>(k0\r(q)  - r ’((?’)]),  the 
tangential  electric  field  can  be  written  in  terms  of  auxiliary  vector  and  scalar  potentials  depending  on 
tangential  surface  current  (Jt)  and  surface  charge  (p)  densities. 
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Esc.t  = m A.t-V( t> ./  = /con  J(/i ' (F'  )G* (F, r')dl'-t.V  0J p (r')G*  (r,r')dl'  (1) 

Then  by  applying  equation  of  charge  conservation,  i.e.,  V^I  — idp,  and  PEC  boundary  condition  to  the 
above  equation  and  after  manipulations  the  following  integral  equation  is  obtained: 


-E'"c (j)P(s)  = /(0H  |j,(5,)cos(^(i)-^(5,))P(s)P(5,)G0/'(5,5')rfe'-^— J—  G*(s,s')ds'  (2) 

where  $(s)=r(s)/cos(i(s)),  E ,(s)  is  the  angle  between  the  contour  normal  and  the  x-direction,  and  y (s)  is 
the  angle  between  the  contour  normal  and  the  radial  direction.  Besides  the  closed  contour  C is 
characterized  by  parametric  equations  x=x(s),  y=y(s),  0<s<2n  and  the  current  density  function  J,  is 
assumed  zero  on  the  part  of  C complementary  to  M.  The  above  IE  is  defined  on  M.  On  the  slot  part  the 
current  density  is  zero.  The  spectral  versions  of  these  two  equations  constitute  a DSE.  However  a special 
care  is  needed  for  the  Green’s  functions  to  find  a regularized  solution.  Therefore  we  define  the  following 
new  auxiliary  regular  functions, 

H(s,s')  = (s’)|)  - //^(2^|sin(5 -s')  / 2|)  (3) 

G(s, s')  = cos($  (s)  - \ {s'))  p {s)  p (s')//(1)(£|F(s)-F(s’)|)  - p 2{s)  H 0)(2fci|sin(s -s')  / 2|) 

(4) 

These  functions  H(s,s  ’)  and  G (s,s  ’)  have  continuous  first  derivatives  and  second  derivatives  with  respect 
to  s and  s’  and  belong  to  L2.  For  regularization  these  functions  are  expanded  into  double  FS  and  the 
corresponding  expansion  coefficients  rapidly  decay  like  0(l/(n\\m\)32'z).  Then  by  using  (3)  and  (4),  DSE 
is  written  in  certain  canonical  form  and  solved  by  RHP  technique.  The  main  problem  in  this  process  is  an 
efficient  computation  of  the  double  FS  coefficients.  We  perform  this  by  using  FFT  algorithms  and  exploit 
the  fact  that  the  second  terms  in  (3)  and  (4)  are  expanded  into  FS  analytically.  All  this  provides  us  to 
handle  larger  problems  than  can  be  treated  with  earlier  methods. 

NUMERICAL  RESULTS:  Firstly  the  convergent  nature  of  our  algorithm  is  shown  in  Figure  2.  The 
relative  error  in  the  computed  surface  current  is  max\xnN,r  ‘ - xnNlr\/max\x„N"\,  that  decreases  quite  rapidly 
with  larger  Nlr.  It  can  be  said  that  0.1%  accuracy  can  be  obtained  with  Nlr=2ka  + 15.  Figure  3 presents  a 
comparison  of  the  radiation  patterns  obtained  by  the  presented  method  and  Physical  Optics  solution.  It 
shows  that  excluding  the  back  lobes  the  two  methods  show  quite  similar  patterns  but  still  some 
disappearances  are  observed  in  the  penumbra  region. 

CONCLUSION:  Two-dimensional  parabolic  reflector  antenna  illuminated  with  directive  feed  is  solved 
by  the  analytical  regularization  in  the  H-polarization  case.  Presented  numerical  results  have  uniform  3- 
digit  accuracy  although  the  algorithm  can  generate  them  with  machine  precision.  Application  of  the  FFT 
codes  provides  us  to  solve  larger  problems,  and  our  data  can  be  used  as  a reference  to  validate  purely 
numerical  or  asymptotic  methods. 
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Figure  1:  Geometry  of  2D  reflector  antenna  system 


Figure  2:  Truncation  Error  in  the 
current  coefficients.  Feed  directivity 
factor  is  &6=2.6 


Normalized  Radiation  Pattern 


Figure  3:  Comparison  of  radiation  patterns 
for  MAR  and  PO  solutions.  D=  20X. 

Feed  directivity  factor  is  kb- 1 .5. 
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